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1  SUMMARY 


This  grant  focused  on  investigating  effects  of  the  atmosphere  on  eleetromagnetie  signal 
propagation  at  W-band  (81-86  GHz)  and  V-band  (71-76  GHz).  We  compared  modeling 
approaches  that  included  dry  air  (gases),  water  vapor,  fog,  elouds,  and  rain.  We  developed 
laboratory  test  faeilities  to  eharaeterize  antenna  performance  at  these  frequencies.  We  developed 
a  slant-path  experiment  to  measure  the  effects  of  weather  and  assess  model  aceuracy.  The  goal  is 
that  through  this  experiment,  a  reliable  attenuation  model  ean  be  constructed  (or  verified)  to 
indicate  the  amount  of  attenuation  a  specifie  terrestrial  link  would  have  at  these  frequencies.  This 
could  then  be  applied  to  satellite  communieation  link  modeling. 

2  INTRODUCTION 

The  W-band  (81-86  GHz)  and  V-band  (71  -  76  GHz)  frequency  bands  have  potential  to 
provide  very  high  bandwidth  (i.e.,  gigabit  per  second  (Gbps)  data  rates),  point-to-point 
eommunieations  between  a  satellite  and  a  ground  station  [1-4].  However,  one  of  the  most  eritical 
problems  in  achieving  this  type  of  eommunication  link  is  the  uncertainty  of  signal  attenuation, 
depolarization,  and  phase  distortion  that  results  from  atmospherie  absorption,  seintillation,  and 
meteorologieal  effeets  (e.g.,  rain- fade,  atmospheric  and  molecular  absorption,  and  scattering 
from  hydrometers).  Rain  parameters  such  as  drop-size,  drop-shape  and  canting  angle  can  impact 
the  propagation  of  waves  at  sueh  extremely  high  frequeneies.  As  rainfall  affeets  the  propagation 
of  the  wave,  rainfall  itself  is  also  affeeted  by  several  other  faetors,  such  as  wind  velocity  and 
temperature  [5-7].  During  a  turbulent  condition,  a  vertical  wind  velocity  component  is  present, 
and  the  raindrop  density  in  the  microwave  beam  does  not  change,  but  the  rain-rate  measured  at 
the  ground  ehanges  considerably  [8-12],  whieh  causes  uncertainty  in  the  eorrelation  of 
attenuation  and  rain-rate  measurements.  Furthermore,  the  overall  dieleetric  eonstant  of  rainfall 
water  is  altered  with  temperature,  thus  ehanging  the  values  of  the  attenuation  coeffieient. 

There  are  several  approaehes  that  seientists  have  used  in  the  past  to  estimate  rain 
attenuation  based  on  the  use  of  a  large  number  of  attenuation  variations  at  different  frequeneies, 
loeations,  and  path  geometries.  Additionally,  there  are  several  predictive  models  that  are  based 
on  theoretical  evaluation  and  extrapolations  to  higher  frequencies  using  lower  frequency 
observations.  Our  goal  in  this  work  was  to  devise  an  experiment  to  colleet  attenuation  and 
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depolarization  data  directly  at  the  WA^-bands  so  high-fidelity  propagation  models  can  be 
developed  and  statistically  validated.  Our  proposed  terrestrial  link  experiment  (simply  referred  to 
as  W/V-band  Terrestrial  Link  Experiment,  or  WTLE)  would  enable  the  measurement  of 
attenuation,  rain-rate,  and  rain-drop  size  distribution  using  the  receiver  and  weather  sensors.  The 
two  main  parameters  of  interests  that  are  expected  to  be  measured  using  this  terrestrial  link  are: 
(1)  rain  attenuation  and  (2)  depolarization. 

Attenuation  is  usually  modeled  by  the  power  law  relationship  given  by: 

A  =  aRp,  (1) 

where  A  is  the  attenuation  due  to  rain,  in  dB/km,  R  is  the  rain  rate,  in  mm/hr,  and  a  and  (3  are 
constants  that  depend  on  parameters  such  as  frequency,  rain,  temperature,  and  rain  drop  size 
distribution  [13].  Depolarization  of  a  propagating  electromagnetic  wave  is  usually  due  to  the 
presence  of  ice  crystals  at  high  altitudes  (typically  greater  than  2-km)  and  rain  at  lower  altitudes 
(less  than  2-km).  In  general,  depolarization  of  a  propagating  signal  can  be  modeled  in  terms  of 
the  scattered-field  matrix  as: 

[EJ  =  (2) 

where  L  is  the  path  length  through  the  scattering  region,  \Ei\  is  the  incident  field  vector,  [E^]  is 
the  scattered  field  vector,  and  and  are  the  attenuation  and  phase-lag  vector  components, 
respectively.  These  components  are  dependent  on  the  rain-drop  geometry  and  distribution.  Both, 
the  attenuation  and  depolarization  constants  vary  statistically  depending  on  local  weather. 
Therefore,  it  is  important  that  the  WTEE  setup  collects  adequate  data  to  enable  studies  of  rain 
attenuation  over  many  rain  events,  with  a  reliable  knowledge  of  rain-rate  and  rain-drop  size 
distribution. 

We  investigated  the  International  Telecommunications  Union  (Radio)  (ITU-R)  models 
for  rain  [14],  atmospheric  gases  [15],  and  clouds  or  fog  [16]  attenuation.  These  models  were 
developed  by  researchers  over  many  years  and  have  been  verified  using  data  collected  world¬ 
wide.  Eor  comparison,  we  investigated  a  model  developed  by  Euis  da  Silva  Mello  that  includes 
path  attenuation  due  to  rain  attenuation  [17].  Comparisons  at  two  frequencies  (72  GHz  and  82 
GHz)  were  made  and  included  in  this  report.  Simulation  results  are  presented  for  the  proposed 
WTEE  experiment.  To  incorporate  realistic  data  for  the  models  of  the  specific  link,  weather  data 
from  local  weather  stations  were  accessed  via  the  Mesowest  website  [18].  Water  vapor  densities 
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were  also  found  through  an  ITU-R  database  [19]  and  a  related  ITU-R  water  vapor  density  paper 

[20], 

Details  of  the  WTLE  test  are  described  in  this  report.  The  WTLE  tests  are  planned  to  be 
conducted  in  Albuquerque,  New  Mexico,  later  in  2015.  The  transmit  assembly  and  receive 
assembly,  along  with  the  related  equipment,  are  described  in  detail.  A  discussion  of  the  antenna 
test  fixtures  and  other  hardware  fixtures  that  we  developed  to  calibrate  the  antenna  is  presented. 


3  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 
3.1  ITU-R  Model 


3.1,1  Atmospheric  Gases 


One  of  the  factors  affecting  the  propagating  signal  is  the  attenuation  due  to  atmospheric 
gases,  such  as  nitrogen,  oxygen,  and  water  vapor.  The  ITU-R  reference:  “Attenuation  by 
Atmospheric  Gases”  provided  a  model  for  these  atmosphere  constituents  [15].  This  reference 
provided  two  different  methods  to  estimate  the  attenuation  on  terrestrial  and  slant  paths.  The 
method  applied  herein  was  a  simplified  approximation  for  the  frequency  range  of  1-350  GHz. 
The  total  (combined)  specific  attenuation  (i.e.,  dB/km)  due  to  atmospheric  gases  was  estimated 
using  a  simplified  algorithm  based  on  curve-fitting  to  the  line-by-line  calculation  of  attenuation 
by  dry  air  and  by  water  vapor. 

Eor  dry  air  (i.e.,  no  water  vapor),  the  attenuation,  yo  (dB/km),  can  be  estimated  by  the 
following  equation: 


To 


=  {3. 


02  X  10 


-4„3.5 
't 


0.283r, 


3.8 

t 


+  X  10-=  (3) 


(/-118.75)2+2.91r2rj^®  +  j" 


where: 


^4  =  (p(rp,rt, -0.0112, 0.0092, -0.1033, -0.0009),  (4) 

^5  =  (p(rp,rt,  0.2705, -2.7192, -0.3016, -4.1033),  (5) 

^6  =  (p(rp,rt,  0.2445, -5.9191, 0.0422, -8.0719),  (6) 

^7  =  (p{rp,  Tt,  -0.1833, 6.5589,  -0.2402, 6.131),  (7) 

f.  frequency  (GHz), 

Tp  =  Plot/ 1013,  where  ptot  represents  total  air  pressure,  (8) 

rt  =  288/(273+0,  (9) 
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p;  pressure  (hPa),  and 
t:  temperature  (°C). 

For  water  vapor,  the  attenuation,  Yw  (dB/km),  can  be  estimated  by  the  following  equation: 


Vw 


3.98??iexp[2.23(l-rt)]  11.96?7i  exp[0.7(l-rt)] 

(/-22.235)2+9.42p^  «v/>  J  T”  (^_i83,3i)2  +  m4^2 


O.OSlpi  exp[6.44(l-rt)]  ^ 
(/-321.226)2+6.29?7^ 


3.66?7iexp[l.6(l-rt)]  25.37?7i  exp[l.09(l-rt)]  17.4?7i  exp[l.46(l-rt)] 

(/-325.153)2+9.22?7^  (/-380)2  (/-448)2 


8«.6,.exp[0.17(l-r^  f,(/,557)  +  .xplo.-ld-r,)]  752) 

v/ 


(/-557)2 


8.3328x10^772  exp[0.99(l-rt)] 

(/-1780)2 


g(^f,nQQi)^prl^p  x 


where: 


r/i  =  0.955rprt°-^^  +  0.006p, 
1/2  =  0.735rprt°-^  +  QmS3r^'^p, 

gif.fi)  -  1  +  (^)', 


(10) 

(11) 

(12) 

(13) 


P 


(14) 


and: 

p:  water  vapor  density  (g/m  ), 

pi:  water  vapor  density  at  altitude  hi  of  the  corresponding  station,  and 
hi:  altitude  of  the  receiving  station. 


These  expressions  can  be  combined  to  estimate  the  total  path  attenuation  due  to  dry  air 
and  water  vapor.  For  WTLE,  propagation  path  is  inclined  at  an  angle  of  about  4°.  Therefore,  the 
procedure  presented  by  the  ITU-R  to  calculate  the  attenuation  of  an  inclined  path  having  an 


elevation  angle  between  0°  and  5°  was  used.  Thus,  the  total  attenuation  can  be  estimated  by: 


A  =  ToV^ 


-Til 

Ao 

-tl2/ 

^Re+h^F(xPe 

^Re+h2F{x2)e  t^o 

coscp^ 

COS(p2 

_  -^2/ 

+h2F(x2')e 

+  yw4^ 


COS(p^ 


COS(p2 


(15) 
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where: 


and 


F(x)  =  - - — 

^  ^  0.661x+0.339Vx2  +  5.5i’ 

(P2  =  arccos(^^-^cos(pi  j, 


Xj  =  tancpi 

x[  —  tancpi 


Re-^hi 


hn 


h\A 


for  i  =  1,2, 


for  i  =  1,2 


(16) 

(17) 

(18) 
(19) 


A:  total  attenuation  (dB), 

Re',  effeetive  earth  radius  (8,500  km),  and 
q)i:  elevation  angle  at  altitude  hi. 

Frequeney,  temperature,  and  air  pressure  are  faetors  required  to  ealeulate  the  attenuation  from 
the  ITU-R  model  equations.  The  altitudes  at  the  referenee  points  are  also  required 


3.1,2  Clouds  and  Fog 

Another  faetor  affeeting  the  propagating  signal  is  attenuation  due  to  elouds  and  fog, 
which  is  distinct  from  “water  vapor”  modeled  in  the  previous  section.  To  model  the  effeet  of 
eloud  and  fog,  the  ITU-R  referenee:  “Attenuation  Due  to  Clouds  and  Fog”  was  used  [16]. 
Aeeording  to  the  ITU-R  paper,  if  the  clouds  and  fog  consist  of  small  droplets  less  than  0.01 -cm, 
and  if  the  frequency  of  operation  is  less  than  200  GHz,  the  Rayleigh  approximation  can  be  used 
to  estimate  the  specific  attenuation  within  a  cloud  or  fog.  Thus,  the  specific  attenuation  (dB/km) 
is  given  by: 

Yc  =  KiM,  (20) 

where: 

Yc:  specific  attenuation  within  the  cloud  (dB/km), 

Ki.  specific  attenuation  coefficient  ((dB/km)/(g/m  )),  and 
M:  liquid  water  density  in  the  cloud  or  fog  (g/m  ). 

To  find  the  specific  attenuation  coefficient,  Ki,  a  mathematical  model  based  on  the  Rayleigh 
scattering  was  used.  This  model  can  be  used  to  calculate  the  value  of  Ki  for  frequencies  as  high 
as  1,000  GHz.  The  equation  for  Ki  is  given  by: 
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where: 


_  0.819/ 

'  f'hi+pZ)’ 

(21) 

2+e' 

V= 

(22) 

""rn  /Go-£i)  I  /(£i-£2) 

fvH  +  if/fpf]  /s[l  +  (///s)1  ’ 

(23) 

(24) 

So  =  77.66+  103.3  (0-  1), 

(25) 

Si  =  0.0671  So, 

(26) 

S2  =  3.52, 

(27) 

0  =  300  /  T, 

(28) 

/p  =  20.20  -146  (0  -  1)  +  316  (0  -  1)^  (GHz), 

(29) 

/s  =  39.8/p  (GHz), 

(30) 

T :  temperature  in  Kelvin;  and 
/:  frequeney  in  GHz, 

Temperature  and  the  liquid  water  density  are  two  faetors  required  by  the  model. 


3.1,3  Rain 


A  third  faetor  affecting  the  propagating  signal  is  the  attenuation  due  to  rain.  To  study  the 
effect  of  that  factor,  the  ITU-R  “Specific  Attenuation  Model  for  Rain  for  use  in  Prediction 
Methods”  was  used  [14].  According  to  the  ITU-R  reference,  the  specific  attenuation  (dB/km)  due 
to  the  rain-rate  can  be  calculated  using  the  following: 

YR^kR^  (31) 

where: 


,  _  [kM+kv+(kif-kv)cos^9  cos  2  t] 

K  —  - , 

2  ’ 

_  [kfjafj+kyay+(kfjafj-kyaY)cos^9  cos  2  t] 

~  2k 

Yr:  specific  attenuation  due  to  rain( dB/km), 

R:  rain  rate  (mm/h), 
k:  constant  coefficient, 
a:  constant  coefficient. 


(32) 

(33) 
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0:  path  elevation  angle, 

x:  polarization  tilt  angle  relative  to  the  horizontal  (x  =  45°  for  eireular 
polarization), 

kn  and  Oh  are  eonstant  eoeffieients  for  horizontal  polarization,  and 
k^,  and  Oy  are  eonstant  eoeffieients  for  vertieal  polarization. 

Values  for  kn,  an,  K  and  Uyfor  different  frequeneies  were  provided  in  the  referenee.  An  estimate 
of  the  rain-rate  was  required  to  ealeulate  the  attenuation. 

3.2  Silva  Mello  Model 

The  path  attenuation  due  to  rain  based  on  the  Silva  Mello  model  will  now  be  explained. 
This  model  used  the  eomplete  rainfall  rate  distribution,  and  the  eoneepts  of  effective  path  length 
and  rainfall  rate.  It  has  been  shown  that  this  method  has  signifieant  improvement  over  the  ITU-R 
method,  beeause  the  ITU-R  “underestimates  the  attenuation  measured  in  terrestrial  links, 
partieularly  for  low  latitude  regions  with  severe  rain  regimes”  [17].  Even  though  Albuquerque  is 
not  eonsidered  to  be  a  region  with  a  severe  rain  regime,  a  eomparison  between  Silva  Mello  and 
ITU-R  was  made. 

For  the  Silva  Mello  model,  the  equation  used  to  ealeulate  the  attenuation  due  to  rain  is 


given  by: 

A  =  (dB),  (34) 

^'^do(Rp) 

where: 

Rgff  =  1.763  j?o-753+o.i97/d  (effective  rain  rate),  (35) 

do  =  (equivalent  eell  diameter),  (36) 

d:  aetual  path  length  (km),  and  (37) 

R:  rain  rate  (mm/h).  (38) 


The  values  of  k  and  a  were  given  by  Olsen  et  al.  [13].  The  parameters  k  and  a  depend  on  the 
frequeney  of  operation  and  the  rain  temperature.  Values  of  k  and  a  for  four  different 
distributions,  whieh  depend  on  the  region  of  operation  and  the  average  drop  size  distribution, 
were  presented  referenee  [13].  These  distributions  were: 
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1)  Laws  and  Parsons  Distribution  (LP):  mainly  used  when  the  rain  rate  is  below  35 
mm/h  in  continental  temperate  rainfall  with  a  mean  drop  size  spectrum; 

2)  Marshall-Palmer  Distribution  (MP):  mainly  used  for  widespread  rain  in  continental 
temperate  climate  (note  that  this  distribution  tends  to  overestimate  the  number  of 
small  drops); 

3)  Thunderstorm  Distribution  (JT):  mainly  used  in  convective  rain  (note  that  this 
distribution  has  not  been  extensively  tested,  but  it  has  been  used  for  some  cases);  and 

4)  Drizzle  Distribution  (JD):  mainly  used  in  regions  having  mostly  small  drops  (note 
that  this  distribution  is  not  generally  used  in  communication  systems  purposes, 
because  it  accounts  mostly  for  small  droplets). 

To  increase  the  accuracy  for  the  Laws  and  Parsons  distribution,  two  regressions  were  performed. 
One  regression  takes  into  account  low  rain  rates  (LPL).  The  other  takes  into  account  high  rain 
rates  (LPH).  Therefore,  in  total,  five  different  distributions  are  available  with  which  to  estimate 
the  values  of  k  and  a.  These  distributions  include:  LPL,  LPH,  MP,  JT,  and  JD. 

The  values  of  k  and  a  were  given  at  specific  frequencies  for  each  distribution.  Since  the 
frequencies  of  interest  (72  GHz  and  82  GHz)  were  not  included  in  the  data  provided,  an 
interpolation  was  performed.  In  the  tables  given  by  Olsen  et  al.  [13],  the  available  data  were 
given  for  specific  rain  temperatures,  which  included  {-10  °C,  0  °C,  and  20  °C}.  Therefore, 
attenuation  calculations  were  performed  only  for  these  rain  temperatures.  The  equations  given  by 
the  Silva  Mello  model  were  used  to  calculate  the  total  attenuation. 

3,3  Weather  Data 

For  each  model,  weather  data  were  required.  For  example,  in  order  to  calculate  the 
attenuation  due  to  atmospheric  gases,  the  temperature,  pressure,  and  water  vapor  density  must  be 
obtained.  Most  of  the  weather-related  inputs  were  obtained  from  the  Mesowest  website  [18]. 
However,  the  water  vapor  density  was  obtained  through  the  use  of  data  provided  by  the  ITU-R 
website  [19]. 

The  Mesowest  website  is  a  database  of  weather  stations  that  allows  users  to  download  the 
weather  data  from  any  of  the  weather  stations  provided.  From  this  website,  the  time  and  dates  for 
the  temperature,  pressure,  and  rain  for  each  of  the  desired  weather  stations  were  obtained. 
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To  calculate  the  attenuation  for  atmospherie  gases,  the  water  vapor  density  was  required. 
This  data  was  not  provided  by  the  Mesowest  website.  The  ITU-R  website  provided  a  database  of 
surfaee  water  vapor  density  values  while  an  assoeiated  paper  [20]  described  how  to  estimate  the 
water  vapor  density  at  any  loeation  on  earth  and  at  any  altitude.  For  the  ITU  model,  the  only 
required  weather  related  input  for  rain  attenuation  was  rain-rate,  whieh  was  provided  by  the 
Mesowest  website.  For  the  Silva  model,  the  required  inputs  were  the  rain-rate  and  the  rain 
temperature.  It  was  determined  that  the  rain  temperature  can  be  estimated  by  the  air  temperature, 
which  was  provided  by  the  Mesowest  website.  The  ITU  model  for  the  attenuation  due  to  clouds 
or  fog  required  the  cloud  thickness  or  the  liquid  water  density.  At  this  point,  a  method  for 
obtaining  this  information  has  not  been  determined.  The  data  was  simply  estimated  based  on  the 
temperature  and  humidity,  with  the  maximum  being  0.5  g/m  for  thick  clouds  and  the  minimum 
being  0.05  g/m  for  medium  elouds. 

To  determine  water  vapor  density,  a  water  vapor  density  paper  from  ITU-R  [20]  was 
used.  The  paper  detailed  the  steps  to  download  the  surface  water  vapor  density  values  on  a  grid 
of  latitude  and  longitude  values  for  speeific  probabilities  of  exeeedanee.  Next,  the  paper 
described  how  to  use  the  data  provided  to  bi-linearly  extrapolate  the  values  for  any  point  on  earth 
and  at  any  altitude.  Using  the  approaeh  mentioned  in  the  paper,  a  eode  was  developed  to 
caleulate  the  water  vapor  density.  The  eode  finds  the  four  elosest  latitude  and  longitude  (by  GPS) 
points  to  the  desired  points.  Onee  it  finds  these  global  positioning  system  (GPS)  points,  it  finds 
the  four  surface  water  vapor  density  values  corresponding  to  the  GPS  points.  It  then  finds  the 
four  vertieal  scaling  factors  corresponding  to  the  GPS  points.  Then,  it  vertieally  scales  the 
surfaee  water  vapor  density  values  based  on  the  altitudes  of  the  four  corresponding  GPS  points, 
the  altitude  of  the  desired  GPS  point,  and  the  four  vertical  scaling  factors.  These  new  four  values 
beeome  the  water  vapor  density  values  at  the  appropriate  altitudes.  Onee  the  vertieally  scaled 
water  vapor  density  values  are  obtained  for  the  four  eorresponding  GPS  points,  the  vertically 
scaled  water  vapor  density,  or  simply  the  water  vapor  density,  was  caleulated  by  bi-linearly 
extrapolating  from  the  four  vertieally  scaled  water  vapor  density  values. 
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3,4  WTLE  Experiment  Description 

The  goal  of  WTLE  is  to  build  an  experiment  test  range  for  72  GHz  and  82  GHz.  The 
objective  is  to  calculate  signal  attenuation  and  correlate  to  atmospheric  factors  and  conditions. 
The  transmitter  would  be  placed  at  Sandia  Peak,  near  Albuquerque,  New  Mexico.  The  receiver 
would  be  placed  at  COSMIAC  (Configurable  Space  Microsystems  Innovations  and  Applications 
Center)  located  in  Albuquerque,  New  Mexico.  This  provides  an  inclined  path  (see  Figure  1) 
where  the  receiver  is  at  an  altitude  of  1.619-km  (5,260  ft)  (above  sea  level)  and  the  transmitter  is 
at  an  altitude  of  3.255-km  (10,670  ft).  The  projected  distance  between  the  transmitter  and  the 
receiver  is  approximately  22.5-km,  and  the  actual  path  length  is  approximately  22.56-km.  GPS 
coordinates  for  the  sites  are  given  in  Table  1.  The  elevation  angle  is: 

=  arctan  ( — — j  ~0.726  rad  ~  4.16  .  (39) 


transmitter 


sea  level 


Figure  1,  Geometry  of  WTLE  Path 


Table  1,  Coordinates  for  WTLE  Receiver  and  Transmitter  Sites 


Receiver  Site 

Transmitter  Site 

Latitude 

35°36.26’N 

35°12.7’N 

Longitude 

106°37.1’  W 

106°27’  W 

Altitude  (ft) 

5,260 

10,670 
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Weather  data  from  three  different  weather  stations  located  between  the  transmitter  and 
the  receiver  along  the  path  would  be  used  for  the  experiment  and  for  simulations.  Figure  2  shows 
the  relative  locations  of  the  primary  weather  stations  along  the  path.  Also  shown  are  the  locations 
of  secondary  weather  station  that  can  provide  additional  data  for  interpolation.  The  KABQ 
weather  station  is  located  near  the  COSMIAC  building.  The  C8328  weather  station  is  located 
near  the  intersection  of  the  1-40  interstate  and  Louisiana  Boulevard.  The  AR227  weather  station 
is  located  near  the  intersection  of  Tramway  Boulevard  and  Tramway  Terrace  Street.  The  distance 
between  KABQ  and  C8328  is  estimated  to  be  6.26-km.  The  distance  between  C8328  and  AR227 
is  estimated  to  be  10.14-km.  The  distance  between  AR227  and  Sandia  Peak  is  estimated  to  be 
6.2-km. 


Figure  2,  Location  of  Weather  Stations 
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3,4,1  Transmitter  Configuration 

The  WTLE  transmitter  assembly  was  acquired  from  Quinstar  Technology,  Inc.  A  block 
diagram  is  provided  in  Figure  3.  The  local  oscillator  (LO)  (10  MHz)  signal  is  split  via  a  divider 
and  routed  to  a  phase-locked  dielectric  resonator  oscillator  (PLDRO)  to  generate  the  intermediate 
frequency  (IF)  (12  GHz)  signal  for  both  the  W  and  V  channels.  The  intermediate  frequency 
signal  is  then  up-converted  to  the  carrier  frequency  and  appropriately  filtered.  Signals  are  routed 
by  waveguide  (W/G)  to  a  circular  polarizer,  and  then  coupled  to  lens  antennas.  A  W/G  power 
detector  was  included  in  the  system  electronics  package  to  provide  a  measurement  of  power 
being  delivered  to  the  antennas. 


Figure  3,  Transmitter  Block  Diagram 

Quinstar  characterized  each  component  in  the  system  and  provided  a  detailed  summary  of 
the  test  results  for  the  electronics  package.  A  summary  of  performance  is  listed  in  Table  2.  While 
some  specifications  were  not  met,  the  package  should  provide  a  stable  and  reliable  beacon  for 
WTFF.  Figure  4  presents  a  picture  of  the  completed  transmitter  electronics  package  with  lens 
antennas.  Figure  5  provides  a  schematic  diagram  of  the  transmitter  package.  The  transmitter 
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would  be  located  in  a  relatively  harsh  environment  and  requires  a  climate  controlled  external 
enclosure.  Figure  6  shows  the  transmitter  package  inside  an  environmentally  controlled 
enclosure. 


Table  2,  Transmitter  Electronics  Performance 


System  Parameters 

Specification 

V-band  Test  Result 

W-band  Test  Result 

Operating  Frequencies 

72  GHz  (nom.) 

84  GHz  (nom.) 

72.000,002,321 

84.000,002,669 

Antenna  Half-Power 
Beamwidth  (HPBW) 

3°  as  a  goal 

3.6°  (E/H) 

3.2°(E)/3.0°(H  ) 

EIRP 

0  dBm/20  dBm  (min.)(manually 
switchable).  3  dB  difference 
between  two  frequency  bands 

41.1 

(  with  5  dB  Attenuator) 

40.4 

(  with  5  dB  Attenuate 

EIRP  Stability 

+/-  0.05  dB(l  sec) 
+/-0.1dB(24hrs) 

-0.041  dBm/  °C 

-0.073  dBm/°C 

Polarization 

Single  LHCP  or  RHCP 

LHCP 

LHCP 

Axial  Ratio 

<1.5  dB 

0.4 

0.29 

Cross-Pol  Discrimination 

>25  dB 

32 

35 

LO  Phase  Noise 

0.1Hz-55dBc 

1.0  Hz  -85  dBc 
10Hz-115dBc 

100  Hz -135  dBc 

-74dBc/Hz 
-107  dBc/Hz 
-130dBc/Hz 
-155  dBc/Hz 

-74dBc/Hz 
-107  dBc/Hz 
-130  dBc/Hz 
-155  dBc/Hz 

LO  Frequency  Stability 

3.00E-12  (<100sec) 
3.00E-9  (per  day) 
l.OOE-7  (per  year) 

1.8E-11 

5E-10 

1.8E-11 

5E-10 

Air  Temperature  Stability 

<0.1  °C  rms 

N/A 

N/A 

Figure  4,  Picture  of  Transmitter  Assembly 
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Figure  5,  Schematic  of  Transmitter  Assembly  (units  in  inches) 


Figure  6,  Transmitter  Assembly  in  Environmentally  Controlled  Enclosure 
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3.4,2  Receiver  Configuration 

The  receiver  site  for  the  WTLE  experiment  is  more  intricate  and  would  include:  (1)  the 
receiver  (Rx)  antenna,  (2)  back-end  radio  frequency  (RF)  components,  and  (3)  several 
atmospheric  measurement  devices.  The  atmospheric  measurement  devices  would  include:  (1)  a 
radiometer  to  provide  an  estimate  of  the  attenuation,  (2)  a  disdrometer  to  measure  rain  drop  size 
distribution  and  velocity,  and  (3)  a  weather  station.  A  component  layout  of  the  receiver  system  is 
shown  in  Figure  7. 

A  block  diagram  of  the  receiver  electronics  package  is  provided  in  Figure  8.  The  flow  of 
the  signal  is,  in  general,  exactly  the  opposite  of  the  transmitter  package.  The  antennas  detect  the 
beacon  signal,  which  is  bandpass  filtered  to  the  respective  frequencies.  The  two  paths  are  co¬ 
polarized  left  hand  circular  polarization  (FHCP)  on  top  and  right  hand  circular  polarization 
(RHCP)  on  bottom.  The  signals  are  then  mixed  with  a  phase-stable  local  RF  source.  A  72  GHz 
reference  signal  is  generated  in  the  receiver  backend,  and  then  routed  through  an  externally 
controlled  switch. 
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Figure  7,  WTLE  Receiver  System 
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Figure  8,  Receiver  Block  Diagram 


3,5  Transmit  Structure  Displacement  Measurement  and  Fortification 

During  initial  Sandia  Peak  site  visits,  a  concern  was  raised  about  the  stability  of  the 
structure  to  be  used  for  the  WTLE  transmit  antenna  assembly.  The  structure  exhibited  a  very  low 
and  under-damped  resonance  frequency.  Therefore,  the  structure  was  stiffened  with  braces  that 
were  secured  to  the  roof  of  the  building  as  shown  in  Figure  9.  Weather  instruments  were  attached 
to  the  reinforced  structure  as  shown  in  Figure  10.  Instruments  included  a  barometric  pressure 
sensor,  rain  gauge  (i.e.,  tipping  bucket),  temperature  and  humidity  probes,  and  wind  monitor. 


Figure  9,  Fortification  of  WTLE  Tx  Antenna  Mounting  Structure;  Before  (Left)  and  After 

Brace  Installation  (Left) 
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Figure  10,  Weather  Station  Installed  on  Transmitter  Structure 


Angular  displacement  measurements  of  the  transmit  (Tx)  and  Rx  anteimas  are  required  to 
enable  us  to  separate  signal  degradations  resulting  from  bore-sight  pointing  error  and  from 
atmospherie  attenuation.  Angular  rate  gyros  were  used  to  eolleet  this  data.  These  angular  rate 
gyros  and  eompanion  data  reeorder  were  also  attaehed  to  the  transmitter  strueture.  Figure  11 
shows  the  inside  layout  of  that  box. 


Power  Supply 

Campbell  Scientific,  ModehPSlOO 


Data  Recorder 

Campbell  Scientific,  Model:  CR-1000 


Angular  Rate  Gyros 
Donner,  Model:  SDG-500 


Figure  11,  External  Enclosure  Internal  Layout 
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Data  was  collected  using  weather  station  instruments  during  the  time-period  of  25-Mar- 
2015  to  lO-Apr-2015.  This  test  also  allowed  us  to  obtain  an  estimate  of  the  magnitude  of 
pointing  errors  to  be  expected  at  the  transmit  site.  The  data  showed  good  correlation  between  the 
wind  monitor  and  the  gyro  outputs,  especially  when  the  wind  magnitude  was  normalized  by  the 
direction  of  the  wind  load  on  the  structure.  The  data  is  presented  in  Figure  12.  The  magnitude  of 
the  recorded  angular  displacement  was  on  the  order  of  a  tenth  of  a  degree.  Unfortunately,  the 
data  exhibited  several  anomalies.  We  are  concerned  that  the  current  location  may  have  very 
strong  electromagnetic  fields  produced  by  neighboring  transmitters,  which  may  corrupt  our 
instruments.  If  this  is  true,  the  WTLE  transmit  electronics  may  not  perform  adequately  and  a 
different  site  may  have  to  be  selected. 


Time  (seconds) 

Figure  12,  Time  Traces  of  the  Wind  and  Gyro  Measurements 

3,6  Anechoic  Chamber  Preparation 

The  WTLE  lens  antennas  are  to  be  measured  and  characterized  by  students  in  the 
Electrical  and  Computer  Engineering  Department  at  the  University  of  New  Mexico  (UNM/ECE). 
Preparation  of  the  UNM/ECE  anechoic  chamber  for  high  frequency  measurements  included  the 
design  of  a  transmit  antenna  mount,  shown  in  Figure  13.  The  mount  provides  a  means  to  ensure 
that  the  output  aperture  of  the  device  under  test  is  centered  directly  over  the  rotating  positioner.  It 
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allows  for  adjustment  in  all  linear  directions.  The  antenna  mount  was  fabricated  out  of  plastic. 
Additionally,  the  hardware  that  secures  the  mount  was  also  fabricated  from  plastic  to  reduce 
interference  with  the  device  under  measurement. 


Figure  13,  Solid  Model  of  Anechoic  Chamber  Tx  Antenna  Mount  Assembly 

Another  important  consideration  for  antenna  measurements  is  the  orthogonality  of  the  Tx 
and  Rx  antennas.  The  antenna  mount  allows  for  fine-tune  adjustment  of  the  y-axis  rotation,  as 
shown  in  Figure  14.  Rotation  in  the  z-axis  is  controlled  by  a  Lindgren  Model  5000  rotation 
positioner.  When  the  antenna  is  at  the  desired  height  in  the  z-axis,  the  nut  on  the  bottom  post  of 
the  yellow  disc  can  be  tightened.  The  nut  on  the  post  provides  a  pivot-point  for  the  yellow  disc  to 
be  rotated  along  the  curved  slot  until  the  antenna  is  perpendicular,  at  which  point  the  upper  post 
is  then  secured. 
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Figure  14,  Tx  Antenna  Mount,  Adjustment  to  Ensure  Tx  and  Rx  Antenna  Orthogonality 

The  antenna  mount  has  a  feature  to  allow  multiple  slices  of  the  antenna  to  be  measured, 
which  allows  for  3-D  antenna  pattern  measurement.  Both  the  yellow  discs  and  a  central  mount 
located  between  the  yellow  discs  have  holes  spaced  at  30°  increments.  Nuts  permanently  located 
on  the  back  side  of  the  yellow  discs  can  be  used  to  secure  a  bolt  from  the  central  mount  to  the 
yellow  disc.  Figure  15,  shows  how  the  holes  in  the  central  disc  (shown  in  transparent  gray)  line 
up  with  the  holes  on  the  yellow  discs. 


Figure  15,  Tx  Antenna  Mount,  Adjustment  for  3-D  Antenna  Measurements 
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A  final  consideration  for  antenna  measurements  is  the  development  of  an  alignment 
system  to  ensure  the  Tx  and  Rx  antennas  are  perpendicular  and  vertically  aligned.  The  system  is 
shown  in  Figure  16.  We  use  a  small  laser  to  accurately  position  the  device  under  test.  Once 
aligned,  the  laser  alignment  system  can  be  removed  from  the  antenna  mount. 


Figure  16,  Laser  Alignment  Setup 


4  RESULTS  AND  DISCUSSION 

This  section  provides  preliminary  simulation  results  for  the  ITU-R  and  Silva  Mello 
model.  Analysis  of  the  models  and  simulations  continues.  However,  some  qualitative 
observations  can  be  made. 

4.1  ITU-R  Models 

The  simulation  results  for  the  ITU-R  models  are  shown  for  72  GHz  and  82  GHz  in  Figure 
17  and  Figure  18,  respectively.  The  water  vapor  density  was  assumed  to  be  7.442  g/m  for  a  1% 
probability  exceedance  during  the  month  of  April  in  the  city  of  Albuquerque,  NM.  In  other 
words,  99%  of  the  month  of  April,  the  water  vapor  density  would  be  7.442  g/m^  or  less.  The 
predicted  temperature,  pressure,  precipitation,  cloud  thickness,  and  attenuation  were  computed 
for  the  three  weather  station  location  previously  described.  The  results  indicate  that  for  the  rain 
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and  cloud  weather  factors,  the  attenuation  at  82  GHz  is  higher  than  the  attenuation  at  72  GHz,  as 
would  be  expected,  sinee  attenuation  due  to  weather  inereases  with  frequency.  As  for  the 
atmospheric  gas  factor,  the  attenuation  at  82  GHz  is  lower  than  the  attenuation  at  72  GHz,  as 
would  be  expected,  since  82  GHz  is  farther  from  the  oxygen  absorption  complex.  The  difference 
between  the  two  frequeneies  decreases  with  a  decrease  in  temperature.  Thus,  the  models  appear 
to  be  working  as  expeeted. 
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Temperature  at  each  station. 
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Attenuation  over  the  entire  path  due  to  precipitation. 


Attenuation  over  the  entire  path  due  to  atmospheric  gases  (H20,  02) 
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Attenuation  at  each  station  from  the  clouds  (liquid  water  density). 
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Attenuation  over  the  entire  path  due  to  the  clouds. 
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Figure  17.  ITU-R  Models  Results  for/=  72  GHz 


Temperature  at  each  station.  Precipitation  the  last  hour  at  each  station.  Attenuation  over  the  entire  path  due  to  precipitation. 
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4.2  Silva  Mello  Model 


Using  the  framework  diseussed  earlier,  representative  variations  of  temperature  and  rain- 
rate  as  a  funetion  of  time  were  simulated  and  are  presented  in  Figure  19.  Temperature  and 
preeipitation  were  eomputed  for  the  three  ground  station  loeations  planned  for  the  WTLE 
experiment.  The  simulations  depict  a  rain  event  occurring  in  the  afternoon.  Station  1  is  clearly 
expected  to  experience  the  greatest  amount  of  precipitation. 


Temperature  art  each  station. 


12-04-2014 


Precipitation  the  last  hour  at  each  station 


12-04-2014 


Figure  19,  Temperature  and  Rain-Rate  as  a  Function  of  Time 

From  the  assumed  (representative)  temperature  and  precipitation  profiles,  the  total 
attenuation  for  the  WTLE  link  was  simulated.  Each  subsequent  simulation  compares  the 
estimated  total  attenuation  for  all  5  distributions  discussed  previously  (i.e..  Laws  and  Parsons 
(LPL  &  LPH),  Marshall-Palmer  (MP),  Thunderstorm  (JT),  and  Drizzle  (JD)).  Eor  each 
frequency,  simulations  were  conducted  assuming  three  different  rain  temperatures  {0  °C,  -10  °C, 
and  20  °C}.  Eigure  20  shows  the  results  for  72  GHz  at  a  rain  temperature  of  0  °C.  Eigure  21 
corresponds  to  72  GHz  and  rain  temperature  of  -10  °C.  Eigure  22  corresponds  72  GHz  and  rain 
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temperature  of  20  °C.  Figure  23,  Figure  24,  and  Figure  25  present  the  simulation  results  at  82 
GHz  for  {0  °C,  -10  °C,  and  20  °C},  respectively.  As  expected,  there  is  significant  fluctuation  of 
the  total  path  attenuation  during  the  rain  event.  In  all  cases,  the  Marshall-Palmer  distribution 
resulted  in  the  greatest  attenuation  prediction.  The  simulations  indicate  that  total  attenuation  can 
exceed  20  dB  for  the  given  rain-rate  at  72  GHz.  The  simulations  also  suggest  that  total 
attenuation  can  approach  30  dB  at  the  given  rain-rate  at  82  GHz.  Further  quantitative  analysis  of 
the  simulation  results  is  underway.  Qualitatively,  it  appears  that  the  Silva  Mello  model  did 
predict  higher  attenuation  than  the  ITU-R  model  for  both  frequencies. 


Total  Attenuatiori  Due  to  Rain 


Figure  20,  Silva-Model  Results  at/=  72  GHz  and  Rain  Temperature  =  0°C 
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Total  Attenuation  Due  to  Rain 


Figure  21,  Silva-Model  Results  at/=  72  GHz  and  Rain  Temperature  =  -10  °C 


Total  Attenuation  Due  to  Rain 


Figure  22,  Silva-Model  Results  at/=  72  GHz  and  Rain  Temperature  =  20  °C 
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Total  Attenuation  Due  to  Rain 


Figure  23,  Silva-Model  Results  at/=  82  GHz  and  Rain  Temperature  =  0  °C 


Total  Attenuation  Due  to  Rain 


Figure  24,  Silva-Model  Results  at/=  82  GHz  and  Rain  Temperature  =  -10  °C 
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Total  Attenuation  Due  to  Rain 


Figure  25,  Silva-Model  Results  at/=  82  GHz  and  Rain  Temperature  =  20  °C 


5  CONCLUSIONS 

This  work  presented  the  well-known  ITU-R  model  for  atmospheric  propagation 
attenuation  and  a  model  developed  by  Silva  Mello.  Each  model  included  absorption  from  gases, 
water  vapor,  clouds,  and  precipitation.  Simulation  results  were  consistent  with  qualitative 
expectations.  The  models  were  applied  to  the  proposed  WA^-band  Terrestrial  Link  Experiment 
(WTLE)  concept.  Simulation  results  suggest  that  low  values  of  rain-rate  (i.e.,  2-3  mm/hr)  could 
produce  more  than  20  dB  of  total  link  attenuation  for  the  WTLE  propagation  path.  As  expected, 
attenuation  increases  with  frequency. 

The  WTLE  concept  to  experimentally  validate  and  verify  our  models  was  described. 
WTLE  consists  of  a  tonal  transmitter  and  companion  receiver  separated  by  22.5  km  and 
including  a  4.5°  inclination.  WTLE  would  collect  meteorology  data  along  the  propagation  path 
and  attempt  to  correlate  signal  attenuation  to  measured  weather  parameters. 

Test  facilities  to  measure  antenna  performance  have  been  established  and  are  functional. 
Likewise,  transmitter  and  receiver  systems  have  been  developed  and  are  currently  being  installed 
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and  tested.  Weather  instruments  are  also  being  installed  and  tested.  To  mitigate  experiment 
uneertainty,  motion  (angular  rate)  sensors  have  been  installed  on  the  mounting  structures,  which 
would  enable  de-correlation  of  pointing  error  from  the  test  data. 

6  RECOMMENDATIONS 

We  recommend  that  additional  modeling  approaches  be  investigated,  which  may  provide 
valuable  insight  on  the  uncertainty  at  our  target  frequency  range.  New  modeling  tools  are 
required  to  extrapolate  available  measurements  to  the  actual  propagation  path  parameters.  We 
anticipate  that  actual  data  collection  on  the  WTLE  system  will  begin  in  August  2015.  We 
recommend  that  data  be  collected  under  many  operation  conditions.  We  recommend  that  the  data 
be  analyzed  to  determine  statistical  trends  and  experimental  uncertainty. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


3-D 

Three  Dimensions 

AC 

Alternating  Current 

BPF 

Band  Pass  Pilter 

°C 

Degree  Celsius 

cm 

Centimeter 

COSMIAC 

Configurable  Space  Microsystems  Innovations 
and  Applications  Center 

dB 

Deeibel 

dBc 

Deeibels  Relative  to  the  Carrier 

dBm 

Decibel  (relative  to  milli-Watt) 

DC 

Direct  Current 

deg 

Degree 

ECE 

Electrieal  and  Computer  Engineering  Department 

E/H 

Electrie  to  Magnetie  Pield  Amplitude 

EIRP 

Effeetive  Isotropic  Radiated  Power 

/ 

Prequeney 

ft 

feet 

gW 

Density,  grams  per  cubie  meter 

Gbps 

Gigabit  per  second 

GHz 

Giga-Hertz 

GPS 

Global  Positioning  System 

HPBW 

Half  Power  Beam  Width 

hr 

hour 

Hz 

Hertz 

IP 

Intermediate  Prequeney 

ITU-R 

International  Telecommunications  Union  (Radio) 

JD 

Drizzle  Distribution 

JT 

Thunderstorm  Distribution 

km 

kilometer 

EHCP 

Eeft  Hand  Circular  Polarization 

ENA 

Eow  Noise  Amplifier 

EO 

Eocal  Oscillator 
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LP 

Laws  and  Parsons 

LPH 

Laws  and  Parsons,  High 

LPL 

Laws  and  Parsons,  Low 

MHz 

Mega  Hertz 

mm 

millimeter 

MP 

Marshal  Palmer  Distribution 

m/s 

Meters  per  seeond 

PCB 

Printed  Cireuit  Board 

PLDRO 

Phase-Looked  Dieleotrio  Resonator  Osoillator 

Rec 

Reoeiver 

RF 

Radio  Lrequenoy 

RHCP 

Right  Hand  Ciroular  Polarization 

Rx 

Reoeiver 

TBD 

To  Be  Determined 

Tx 

Transmitter 

UNM 

University  of  New  Mexioo 

USB 

Universal  Serial  Bus 

W/G 

Wave  Guide 

WILE 

W/V-band  Terrestrial  Link  Experiment 
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